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Abstract 

The effect of short wavelength quantum electrodynamic (QED) correction on plasma-wave prop- 
agation is investigated. The effect on plasma oscillations and on electromagnetic waves in an un- 
magnetized as well as a magnetized plasma is investigated. The effects of the short wavelength 
QED corrections are most evident for plasma oscillations and for extraordinary modes. In par- 
ticular, the QED correction allow plasma oscillations to propagate, and the extra-ordinary mode 
looses its stop band. The significance of our results is discussed. 
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I. INTRODUCTION 



The non-linear properties of the quantum vacuum has become increasingly interesting 
(see e.g. Refs. for some recent examples), in particular due to the rapidly 

growing power of present day laser systems 0, 0]. It is expected that already the next 
generation laser systems will reach intensities where quantum electrodynamic (QED) effects, 
such as electron-positron pair creation and elastic photon-photon scattering, will be directly 



observable 
experiments 



10, 



_ These effects may even play an important role in future laser-plasma 
^J, e.g., in laser self focusing [l^ ll^ where laser pulse compression could 
give rise to field strengths close to the critical field strength, Ecru = rn^c^/he ~ 10^^ V/m, 
where m is the electron rest mass, c is the speed of light in vacuum, h is Planck's constant 
divided by 27r, and e is the magnitude of the electron charge. Furthermore, in astrophysical 
environments highly energetic phenomena may give rise to parameter ranges in which QED 
can be influential. One such example is pulsar magnetospheres and magnetars 3l- In 
the latter the magnetic fields can reach values above 10^^ G, thus giving equivalent electric 
fields above the Schwinger limit. Thus, in such highly magnetized environments any plasma 
would be one-dimensional and local gamma-ray sources may be prolific. 

These QED effects arise due to the intense electromagnetic field interacting with the 
quantum vacuum. By contrast, in this paper we have chosen to look for linear QED effects 
in plasma wave propagation when the intensity is low. Instead, the frequency is high such 
that short wavelength corrections may be important. The study is carried out for both 
plasma oscillations and for electromagnetic waves in an unmagnetized as well as a magnetized 
plasma. For high plasma densities or high frequency waves, it is found that the short 
wavelength QED corrections can give appreciable corrections to plasma oscillations and the 
propagation of extraordinary electromagnetic modes. We discuss possible applications of 
the results. 



II. BASIC EQUATIONS 

An effective theory for photon-photon scattering can be formulated trough the 
Heisenberg-Euler Lagrangian density Q], which describes a vacuum perturbed with 
a constant electromagnetic field. It is valid for field strengths much lower than the critical 
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field strength and for oscillation length scales much longer than the Compton wavelength. 
By adding derivative correction terms to the Lagrangian density, the effect of rapid field 
variations can be taken into account However, electron-positron pair creation is not 

included in this model, and thus the frequencies must still be much lower than the Compton 



frequency |l2,|2Q| i^e = ttic^ /h, and the dispersive and diffractive effects must be small 21 1. 



The Lagrangian density is 
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where Cq is the classical Lagrangian density, while Che represents the Heisenberg-Euler 
correction due to first order non-linear QED effects, £_d is the derivative correction, F"'^ 
is the electromagnetic field tensor, F"'' = e"'^'^'^F(.d/2 and □ = dad"'. The parameter k = 
2a^h^ / A'bm^c' gives the non-linear coupling, the parameter a = {2/15)ac'^/uj'^ gives the 
dispersive effects of the polarized vacuum, and a = e"^ /Ankceo is the fine structure constant, 
where eo is the free space permittivity. By introducing the four potential such that Fab = 
dgAh — d h An , we obtain the field equations from the Euler-Lagrange equations db [dC/dFab] = 
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where is the four-current and /iq is the free space permeability. 

We further assume that a high frequency low amplitude field is considered, such that 
the non-linear coupling can be neglected compared to the dispersive effects, i.e. the terms 
containing k is negligible compared to the terms containing a. To be more precise, we require 
that eQK\Fab\'^ ^ olo^ i.e. \Fab\lEcr\t ^ ^l^e^ where lo is the frequency with which the 
field typically changes. The corresponding Maxwell equations resulting from the derivative 
corrected field equation then become 
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where p is the total charge density and j is the current density, while the source-free Maxwell 
equations read V ■ B = and 

~~dt' 



V X E 



(2c) 



The fluid continuity and force equations become 
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+ V ■ KVe) = 0, (3) 



dt 
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, V,-V V, = (E + VeXB), (4) 

at J m 

where Ue is the electron density and Ve is the electron fluid velocity. In order to clearly distin- 
guish the QED effects, thermal effects will be neglected throughout this paper except in the 
discussion of plasma oscillations where they will be compared to the QED corrections. Since 
a finite temperature introduces dispersive effects for plasma oscillations, thermal plasma dy- 
namics will compete with the QED corrections (See the Appendix). We note that thermal 
effects may be important in comparison to the QED corrections for other wave modes, but 
this is a topic for future research. 

Since we focus on high frequency phenomena, the ion-motion is omitted. The charge 
density and the current density can then be written as 

p = -e(ne - riio), (5) 
j = -en^Ve, (6) 

where n^o is the constant background ion density. 



III. PLASMA OSCILLATIONS 

Because of the simple classical dispersion relation for plasma oscillations, u"^ = for a 
cold plasma, the dispersive effects due to the QED derivative correction are easy to recognize. 
Below, we neglect thermal effects in our calculations. A complementary discussion where 
thermal effects are included in plasma oscillations is found in the Appendix. Linearizing 
Eqs. pajl and (jH))^© and Fourier decomposing gives us the dispersion relation 

cu^ - cuj = 0, 

where ujp = \/ e^n^j e^m is the plasma frequency. We normalize the parameters according to 

a=^^a, K = uj = —. 

& UJ^ UJp 

The dispersion relation is then written as 

[1 - a {K^ - uj"^)] cu^ - 1 = 0. 
4 



The deviation from the classical result is illustrated in Fig.Q] It is seen that the for small wave 
numbers, the solution follows the classical one (solid line), but the deviation becomes 
stronger as K increases. It is interesting to note that for sufficiently large the phase 
velocity approaches that of the speed of light in vacuum, whereas in the classical case there 
would be no propagation at all. However, within this model the QED corrections must 
remain small such that the condition a^K'^ — o)^) ^ 1 holds. 
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FIG. 1: Dispersion diagram for a plasma oscillations in a cold plasma for different values of a. The 
solid line is the classical dispersion diagram. 



IV. ELECTROMAGNETIC WAVE IN AN UNMAGNETIZED PLASMA 

Linearizing Eqs. ()2b|) . (j^cj, (|H) and © and Fourier analysing, the normalized dispersion 
relation is 

[1 - a [p - [y? -K')-\=^. (7) 

Using the normalized phase velocity v, where v = oj/K = v^p/c, where v^p is the phase 
velocity, we can write the dispersion relation as 

[1 - {atu^ + 1) Lo^] {1 + 2auj'^) u?v^ - auj^= 0, 

which has two distinct solutions. From Fig. |2l it is clearly seen that one of the solutions 
gives the classical limit when a — * 0. For small a;, the other root can approximately be 
written as ~ au'^ to the lowest non-vanishing order. But = tu^/A'^, which implies that 
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~ l/o", and this violates the condition (j{K'^ — uo"^) <^ 1, i.e., the dispersive effects are no 
longer small. For this reason, the non-classical root to the dispersion relation is found to be 
a non-physical one. 

a =0.05 




0.5 1 1.5 2 2.5 



CO 

FIG. 2: The = v^/c? vs. ijd dispersion diagram for an electromagnetic wave in an unmagnetized 
plasma. The dotted line is a non-physical solution to 0. The solid line is the classical dispersion 
diagram. 

The effect of the corrections due to a rapidly varying field is to dislocate the cut-off 
frequency to a slightly lower frequency. The shift in frequency is very small and does not 
become significant until a is much larger than allowed in our model. Thus, we can say that 
within this model, an electromagnetic wave in a unmagnetized plasma is virtually unaffected 
by the short wavelength QED correction. 

V. ELECTROMAGNETIC WAVE, B ||k 

Next we look at wave propagation parallel to an external magnetic field, Bq = -Bqz, where 
the electric field is circularly polarized. 

El = i?i(x ± iy) exp{ikz — iut). 

Linearizing Eqs. ()2b|). (jjcj, (0)) and (jHI) we find the normalized dispersion relation to be 

1 = (C^ ± 0) [1 - CT - u^)] [cj - K^/u}) , 
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where Q = Uc/ujp, uo^ = —eBo/m is the electron gyro frequency and the +(-) sign means a 
right hand circularly polarized wave (left hand circularly polarized wave), also called i?-wave 
(L-wave). The dispersion relation can also be expressed in terms of the normalized phase 
velocity v according to 

[1 - {acu^ + l){iu± n) iu] v^+{l + 2acu^) {to ± 0) cov^ - auj^ [uj±0) = 0. 

In the limit of no magnetic field (fi 0), the dispersion relation simply becomes that 
of an electromagnetic wave propagating in a unmagnetized plasma, which had a non- 
physical root. The corresponding non-physical root has been disregarded in Fig. Eland Fig. 

H 

The effect of the corrections due to a rapidly varying field is again to dislocate the cut-off 
frequency of the R-wave and the L-wave to a slightly lower frequency. Even though the 
effect is more pronounced for the R-wave than for the L-wave, see Fig. IHl and Fig. |3J it is 
still very small. 
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FIG. 3: The = v'^/c^ vs. u) dispersion diagram for an L-wave. The solid line is the classical 
dispersion diagram. 
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o =0.05,£2 =-1 
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FIG. 4: The v"^ = v^/c^ vs. cj dispersion diagram for an i?-wave. The soHd hne is the classical 
dispersion diagram. 



VI. ELECTROMAGNETIC WAVE, B_Lk (X WAVE) 

We now consider an electromagnetic wave propagating orthogonal to the magnetic field 
of a magnetized plasma, such that the electric field oscillates orthogonal to the external 
magnetic field (X-wave). Linearizing Eqs (j2bj) . (f2c|l . (ji} and (jH)) we find the normalized 
dispersion relation to be 

In the limit of no magnetic field (f2 — 0), the dispersion relation becomes a product of that 
for plasma oscillations and that for an electromagnetic wave in a unmagnetized plasma. The 
corresponding non-physical root of (jH)) has been disregarded in Fig. El 

The solutions of (jH)) are easily identified with the two classical ones (solid lines) in Fig. |31 
As can be seen, the high frequency branch does not significantly deviate from the classical 
one. However, the deviation is more pronounced for the low frequency branch. While it is 
not explicitly shown in Fig. El it should be noted that the QED induced effects depend only 
weakly on the external magnetic field. In Fig. jSl this solution has been plotted for different 
values of the normalized dispersive parameter, a. It is interesting to note that the short 
wavelength correction removes the stop band of the X-wave. 
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FIG. 5: Dispersion diagram for an X-wave for different values of a. The solid line is the classical 
dispersion diagram. 



VII. SUMMARY AND DISCUSSION 

In this paper we have investigated how the short wavelength QED correction affects 
plasma wave propagation in both a unmagnetized and a magnetized plasma. In order to 
concentrate on the effects associated with the QED corrections, we have chosen to consider 
the simple case of a cold plasma. Furthermore, the field amplitude is assumed to be small 
such that the nonlinear effects can be neglected. The dispersive effects due to the short 
wavelength QED correction is found to be small but well pronounced for plasma oscillations. 
Fig. and X-waves, Fig. for sufficiently large wavenumbers, whereas the effect is less 
pronounced for parallel propagating R and L-waves, Fig. EJand Fig. El and electromagnetic 
waves in a unmagnetized plasma. Fig. |21 We know that the condition (t(^^ — lj"^) ^ 1 and 
uj Ue must be satisfied, for our model to be apphcable, with a ^ (1/500) (cj^/c^g). For all 
of the waves considered above, except for i?-waves, we learn from the dispersion diagrams 
that u! > Up always, at least for sufficiently large wavenumbers K. Thus, for sufficiently 
dense plasmas, we will reach parameter ranges where the dispersive QED effects become 
large. The highest plasma densities on earth are found in laser fusion experiments where 
the density can reach values as high as Uq ~ 10^^ cm~^ j^], giving Up ~ O.OlcUg. For such 
plasma densities, the dispersive parameter becomes, a ^ 2 x lO"'^. Detection of these QED 
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corrections must therefore be done in the extreme short wave hmit ^ 1 or with very high 
precision. The phase shift of a laser beam can be measured with extremely high precision 



through interferometry [2J|. Also, by using a experimental setup consisting of a standing 
wave in a plasma, the phase shift at each pass through the plasma can be accumulated, thus 
enhancing the total phase shift to be measured. However, to single out the dispersive QED 
effects from the general dispersive effects of a plasma requires detailed knowledge of the 
plasma parameters for such an experiment to be conclusive. In this respect the conditions 
in laser fusion experiments are not suitable, due to the rapid change of plasma density. For 
lower plasma densities, a fairly high stability of the plasma conditions could compensate for 
the even smaller value of a to make a high precision experiment possible. In an experiment, 
the plasma parameters may either be determined with independent measuring techniques, 
but it might also be possible to extract information about the plasma by measuring the 
phase shift for different frequencies or polarizations. In this case, thermal effects may be of 
importance (see the Appendix), and the inclusion of such is therefore of interest for future 
research. How to extract the phase information and how to construct a scheme to detect 
the QED short wave corrections is nontrivial and requires extensive work. 
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VIII. APPENDIX 

In order to highlight the QED effects we have considered the simple case of a cold plasma 
in the manuscript. However, since the QED effects in most cases are small corrections, 
typically even very low temperatures give a comparatively large contribution to the dis- 
persion relations studied, except for wave modes that have density perturbations that are 
identically zero. To illustrate the relative contribution from QED and thermal effects, we 
reconsider the case of plasma oscillations by adding a pressure term — VP/mrig to the right 
hand side of Eq. (j^. Proceeding as in section HI, we linearize the pressure term according 
to — VP/mng = —v^Vne/uQ (where deviation from isothermal behavior is incorporated in 
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the definition of vf), and tlie dispersion relation is then modified to 

,2 



2a 



UJ 



which (given aiuJ^ jc^ — k^) <^ 1) can be approximated by 
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In the regime where Landau damping is small, kvt <^ ujp, we can omit the therm proportional 
to both QED and thermal corrections and write to first order 



UJ 



UJZ 



2aujl 



k' {v^ + 2aujl) . 



Thus the QED effects give a relative frequency shift 2auj'^/c^ of the plasma frequency, and a 
modification of the group velocity of the Langmuir waves proportional to 2auj'^lv\. We note 
that the group velocity shift is a small parameter for any reasonable plasma temperature. 
This suggests that in order to avoid confusion with small temperature fluctuations in a 
possible detection scheme of the QED effects, it might be preferable to look for a small shift 
in plasma frequency, rather than the group velocity shift. 
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